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A gauge invariant theory for electronic Raman scattering for superconductors 
with anisotropic pairing symmetry is analyzed in detail. It is shown that Raman 
scattering in anisotropic superconductors provides a wealth of polarization-dependent 
information which probes the detailed angular dependence of the superconducting 
ground state order parameter. The Raman spectra shows a unique polarization de- 
pendence for various anisotropic pair-state symmetries which affects the peak posi- 
tion of the spectra and generates symmetry dependent low frequency and temperature 
power-laws which can be used to uniquely identify the magnitude and symmetry of 
the energy gap. In particular, we calculate the collective modes and the subsequent 
symmetry-dependent Raman spectra for a d x 2_ y 2 superconductor and compare our 
results to the relevant data on the cuprate systems as well as theoretical predictions 
for s-wave, anisotropic s-wave and s + id energy gaps. Favorable agreement is shown 
with the predictions for d x i_ y i pairing and the experimental data on YBa2Cus07, 
Bi2Sr2CaCu208 and T^B^CuOg. 



(Received February 1, 2008) 



Typeset Using REVTEX 



1 



I. INTRODUCTION 



Knowledge of the symmetry of the energy gap in superconductors provides a major 
step towards unraveling the puzzle of superconductivity in unconventional superconduc- 
tors. While the evidence continues to accrue, the identification of the pair state in the 
heavy-fermion and cuprate systems has proven to be somewhat elusive |T|. While recent 
photoemission experiments @ do allow for an angular-dependent determination of the gap 
and Josephson tunneling measurements have probed the phase of the gap ||, by far the 
most abundant information that has led to speculation of a non-BCS ground state has 
been focussed on the presence of power-laws in the low frequency and/or temperature be- 
havior of transport and thermodynamic quantities ||]. However, due to the averaging over 
the entire Fermi surface, it is well known that the power-laws themselves do not uniquely 
identify the ground state symmetry of the order parameter but only can give the topology 
of the nodes of the energy gap along the Fermi surface, e.g., whether the gap vanishes on 
points and/or lines on the Fermi surface. Thus one cannot distinguish between different 
representations of the energy gap which have the same topology. For instance, for the case 
of o?-wave tetragonal superconductors, there are five pure representations which have line 
nodes on the Fermi surface. Further, the energy gap can become smeared due to inelastic 
quasiparticle collisions, making a fully gapped superconductor seem like a one with gap 
nodes. While in principle the latter effect can be minimized by limiting the experiment to 
very low temperatures, the two-particle correlation functions determining the density, spin 
or current responses do not have the freedom to probe various portions of the gap around 
the Fermi surface, presenting a fundamental obstacle to uniquely identifying the pair state 
representation for the superconductor. 

However, it is well known that Raman scattering has the ability to measure various 
degrees of freedom by simply rotating the incident and scattered photon polarization ori- 
entations. Formally, Raman fluctuations may be viewed as anisotropic mass fluctuations 
around the Fermi surface which do not obey a conservation law [H, such as, e.g., density 
fluctuations. While the isotropic density fluctuations between unit cells will be screened 
due to their coupling to long-range Coulomb forces, the intracell mass fluctuations can 
be anisotropic around the Fermi surface with no net charge and thus can be unscreened, 
providing a scattering mechanism for incoming photons. 

The intensity of scattered light in a Raman experiment can be written in terms of a 
differential photon scattering cross section d 2 a/duidVL as 
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r S 77 (q,^), 



duodVt oj\ 

S 77 (q, w) = [1 + n(u)] Im x 77 (q, w)- (1) 

7T 

Here r = e 2 /mc 2 is the Thompson radius, cji, co>s is the frequency of the incoming and scat- 
tered photon, respectively, and we have set % = ks = 1- 5 , 77 is the generalized structure func- 
tion, which is related to the imaginary part to the Raman response function x 77 through the 
fluctuation-dissipation theorem, the second part of Eq. (1). Finally n(u) = l/[exp(w/T) — 1] 
is the Bose-Einstein distribution function. The generalization of the usual density-density 
correlation (S) and response (x) functions to the case of Raman scattering can be generated 
by weighting the sums over momenta k with the square of the Raman vertex 7(k). 

For small momentum transfers and incident light energies smaller than the optical band 
gap (non-resonant scattering), the vertex for Raman scattering can be written in terms of 
the curvature of the energy band dispersion e(k), 

TtkH^g^, 0) 

where e SI denote the scattered and incident polarization light vectors, respectively, which 
select elements of the Raman tensor. The symmetry of the underlying crystal can be taken 
into account by expanding 7 in terms of a complete set of crystal harmonics </>£, defined on 
the Fermi surface, i.e., [[]] 

7(k) = E7rC(k), (3) 

L,m 

where the index L represents the L-th order contribution to the vertex which transforms 
according to the m-th irreducible representation of the point group symmetry of the crys- 
tal. The quantum numbers L, m classify the anisotropy of the Raman fluctuations around 
the Fermi surface. The full fc-dependence is thus described by the addition of the basis 
functions (which become progressively more anisotropic for higher L ) with different weights 
7™. While the charge, spin, and current density response probes only a single L channel 
(7(k) = 1, L = 0, for the charge and spin density, while 7(k) = k, L = 1, for the current 
density), in principle all even L channels can contribute to the Raman vertex for bands 
which are non-parabolic. Thus by choosing the polarization light vectors accordingly, one 
can select different L, m channels which allow for different projections onto the Fermi sur- 
face. Therefore, electronic Raman scattering can provide detailed information which allows 
insight to the angular dependence of the pair-state symmetry. 
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In accord with this fact, the experimental results on the cuprate systems reveal a wealth 
of polarization-dependent information that provides detailed evidence for determining the 
actual symmetry of the gap 0. The existing body of data on the cuprate systems || reveal 
five main points: 1) in contrast to conventional superconductors such as NbsSn, no clear 
well-defined gap is seen for any polarization orientation even at the lowest temperatures 
measured (T/T c = 0.03) 0, 2) the peak of the spectrum lies roughly at 30% higher frequency 
shifts for the polarization orientation which selects Bi 9 symmetry compared to all other 
symmetries [10-13], 3) there are indications that the temperature dependence of the peak 
in the Bi 9 channel follows more closely to a BCS form than any other symmetry [fL^I , 4) 
the low frequency Raman shifts vary roughly as uj 3 for Bi 9 symmetries and linearly in u for 
the others [9-13], and 5) the ratio of residual scattering in the superconducting state to the 
normal state is smallest for the Bi 9 case compared to all other configurations |T3" |. Such a 
rich spectrum of information should provide a stringent test for various candidates of the 



pairing symmetry states fll4" |. 

The purpose of the present paper is to investigate the polarization dependence of the 
Raman spectra for a superconductor in the weak coupling limit with anisotropic pairing 
symmetry. We calculate the electronic Raman scattering for a tetragonal superconductor 
at finite temperatures and for various polarization orientations in a gauge invariant manner 
and find a rich polarization dependence of the spectra that can be used to uniquely identify 
the energy gap. In particular, we examine the Raman response for s-, d-, s + id-, and 
anisotropic s-wave superconductors and find that a d x 2_ y 2 state agrees surprising well with 
the current information on electronic Raman scattering in YBa 2 Cu3C>7, Bi 2 Sr 2 CaCu2 08 and 
Tl 2 Ba 2 Cu0 6 . 

The plan of the paper is as follows: The ground work for the gauge invariant theory of 
electronic Raman scattering in anisotropic superconductors is reviewed in Section 2 using a 
kinetic equation approach. Section 3 concerns the connection of band structure to the Raman 
vertex and its relation to Fermi surfaces. Section 4 gives our results for the Raman response 
evaluated on a cylindrical Fermi surface for four types of energy gaps, i) s-wave, ii) d-wave, 
iii) s+id-wme, and iv) a special type of anisotropic s-wave. Section 5 presents a comparison 
of the theory to the data on three cuprate systems and contains our conclusions. Appendix 
A deals with a solution of the weak coupling gap equation for the case of anisotropic gaps, 
Appendix B contains details of the theory for the case of (elastic or inelastic) quasiparticle 
scattering in the normal state, while lastly Appendix C is devoted to our calculations for the 
massive collective modes and the subsequent role of vertex corrections using a diagrammatic 
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approach. A brief account of our work has recently appeared 0. 

II. KINETIC THEORY OF ELECTRONIC RAMAN SCATTERING IN 
ANISOTROPIC SUPERCONDUCTORS 

A. Formalism 

In this section we describe a kinetic equation approach for calculating generalized gauge- 
invariant response functions in anisotropic superconductors, with the majority of our atten- 
tion being focused on the electronic Raman response. Other formulations of the calculation 
are possible, and in particular a diagrammatic approach for the gauge invariant Raman re- 
sponse is discussed in detail in Appendix C. Since this approach has been discussed before 
by us [0], we will concentrate on the kinetic equation approach here. 

We consider an anisotropic superconductor in which the electronic states are character- 
ized by a momentum k, an energy (band) dispersion e& = /i + (with \x the Fermi energy), 
a (band) group velocity = V/^, an inverse effective mass tensor (k) = d 2 e k /dkidkj, 
an energy gap A k and excitation energies E k = + |Afc| 2 ] 1//2 . In global thermodynamic 
equilibrium such a system is described by a diagonal equilibrium phase space distribution 
function n k 

n k =<c{c k >=ulf(E k ) + v 2 k [l-f(E k )] = ^[l-^e k ] , (4) 

with 9 k = (l/2E k ) tanh(E k /2T), f(E k ) the Fermi function taken at the Bogoliubov quasi- 
particle energy E k and the usual coherence factors u\ = |[1 + £, k /E k ] and v\ = |[1 — ^ k / E k ] 
for particle-like and hole-like Bogoliubov quasiparticles. In a superconductor there exists 
in addition to the diagonal average n k the off-diagonal average g k 

g k =< c_ k c k >= -A k 9 k . (5) 

The superconducting equilibrium energy gap is then determined from the self consistency 
equation 

A fc = (6) 
p 

with V kp = — \V kp \ the pairing interaction. 

Such a system is assumed to be subject to external perturbation potentials U k ext (q, uj) oc 
exp(iq- r — iut) which are generated in the usual way by an expansion of the Hamiltonian oc 
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(p — eA ext /c) 2 /2m + e$ ext containing q and uj dependent scalar and vector electromagnetic 
potentials $ cxt (q, u) and A ext (q, uj), respectively, to second order in the vector potential 
(including p • A-terms): 

U k ext = e$ cxt + v fe ■ (-^A cxt ) + ^ A\ ■ M^(k) ■ A% 

= e $« rt + Vfc .5ei+7(k)5e 7 , (7) 

where 

8e x = --A cxt , 

c 

Se-y = r |A ! | |A S |, 

7(k) denotes the Raman vertex given by Eq. (2). The last term in Eq. (7) can be interpreted 
to describe the scattering of an incident photon of frequency uj\ represented by the vector 
potential A 1 , into electronic excitations such as particle-hole, Bogoliubov quasiparticle or 
magnon pairs and an outgoing photon (Stokes process) of frequency ujq, = oj\ — uj (vector 
potential A s ) with an associated total momentum transfer q. 

The (linear) response of the distribution function Sn k (q, a;) = n fc (q, uj) —n° k in the absence 
of dissipation is given by the solution of the collisionless (particle-hole-symmetric) kinetic 



equation 1T5L 16 



5n k = —^ — (4> k - X k )5e k - \ k 5el + \ k {u + r] k ) ^— - , (8) 
u)-r] k 2|A fc | 

where t] k = v k - q, 

dnl \A k \\ , e 

h 



and tp k = —df(E k )/dE k . In addition we have defined 



2 A 



k\ 



The cases s — — 1 and s — +1 distinguish the coupling of the response of the pair-correlated 
electron system to phase and amplitude fluctuations of the order parameter, respectively, 
and will be discussed later. The quantity X k is the pair response function introduced by 



Tsuneto |17| which, in the limit q£ <C 1, with £ the coherence length, is given by, 

A fe (q,u;) - -4|A fe | . (11) 
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A more general expression for of valid also for £ 1 < q <C kp is given in Ref. fl8|. The 



total diagonal quasiparticle energy shift 5e k = Se^ + 5e^ may be decomposed into the sum of 
contributions even (+) and odd (-) with respect to the parity operation k — > — k. It differs in 
general from the contribution from the external potential Uj^ through vertex corrections, or 
more physically, through electronic polarization potentials [[19]. This fact may be expressed 
through the following diagonal self-consistency relation, 

6e k = + 2 J2(V q + fk P )Sn p . (12) 
p 

Here, V q = 4-7re 2 /g 2 is the Fourier transform of the long-range Coulomb interaction and f kp 
denotes the short-range Fermi-liquid interaction the consequences of which are, however, 
not considered in what follows since we are interested in the limit of not too large q where 
the polarization correction from the long-range Coulomb interaction dominates the diagonal 
energy change. 

The kinetic equation for the linearized off-diagonal distribution function Sg k (q, to) = 
flffc(q,w) - gl reads 

Sg k + 6 k 5A k + ^^l\ k 5A k - X k SA^ = Ue+ + Vk 5e^} . (13) 

4|Ajfc| \l\ k \ Z\l\ k \ J 

A straight-forward variation of the equilibrium gap equation (6) leads to an off-diagonal 
self-consistency relation 

5A k (q, w) = V k P %(q, £*>), ( 14 ) 
p 

which can now be used to compute the off-diagonal energy shifts , namely the order pa- 
rameter fluctuations 5A k (q,u) and SA\(— q, —uo). They represent the collective oscillations 
necessary to maintain gauge invariance, and must be determined self consistently with the 
off-diagonal kinetic equation |i5| , |T6|| (from which particle-hole asymmetric terms, which are 
typically of the order 0(T/Tp), with Tp the Fermi temperature, have been omitted), 

E V^-rr^fSA; = E V kp X p P ' P . (15) 

In deriving (15) the equilibrium gap equation (6) has been used for simplification. In case 
of particle-hole symmetry, the density, current and Raman fluctuations do not couple to 
the quantity 5A + , which represents the amplitude fluctuations of the order parameter. The 
physical significance of the quantity SA k becomes clear in the (macroscopic) limit ui 2Aq, 
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in which only fluctuations of the phase of the superconducting energy gap determine the 
dynamics of the condensate and <5AjT (q, ui) = \A k \i5(fi(q,uj). 

The important macroscopic observables, namely the density fluctuations 5ni(q,uj), 
charge fluctuations Sn e (q, uj), the Raman fluctuations 5n 7 (q, uj) and the current fluctua- 
tions j(q, uj) are defined as 



8m{q,uj) = 2^1 8n k (q, 



k 



5n e (q,uj) = 2^2e 8n k (q, 



UJ 



UJ) 



k 

<5n 7 (q, uj) = 2 lk8n k (q, uj), (16) 

k 

j(q, uj) = 2 ^ v k [5n k (q, uj) + <p k 5t k (q, uj)]. 

k 

The factors of 2 arise from spin degeneracy. 

We proceed with a solution of Eq. (15). For the time being we would like to restrict 
ourselves to the case where the pairing interaction factorizes as 

V kp = ~V^. (17) 

This ansatz is sufficiently general to allow for an equilibrium gap function A k of arbitrary 
anisotropy in k-space. The maximum of such a gap is denoted A . A k is determined from 
the following form for the equilibrium gap equation (see Appendix A for further details), 

77= E OM-. (18) 

Using Eq. (17), Eq. (15) can be solved immediately to give 

8A~(q,uj) i Y. v K(uj5e+ + r] v 5e-) 

_ 5j5 - j _ = ^( q , u )= ^ . (19) 

The physical significance of the result (19) is that it describes the Goldstone mode for 
superconductors, the Anderson-Bogoliubov or gauge mode, i. e. the massless collective 
mode related to the spontaneously broken gauge symmetry. It is the existence of this 
mode which guarantees gauge invariance of the response theory or, equivalently, charge 
conservation, which we would like to briefly demonstrate now. For this purpose let us write 
the kinetic equation (8) in a form in which the l.h.s is reminiscent of the usual Landau-Silin 
equation |I9| , 



dn° 5A~ 
u5n k - q ■ v k [5n k - ~^$e k ] = -X k [uj5e^ + r) k 5e^] + X k [uj 2 - r)l]— (20) 
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Inserting Eq. (19), one observes that the r.h.s. of Eq. (20) vanishes upon summation over 
momenta k and the l.h.s. of (20) generates the continuity equation for the particle number 
density, 

cj5n 1 -q-j = 0. (21) 

Hence we have demonstrated that accounting properly for the fluctuations of the phase of 
the order parameter leads to the conservation law for the (charge) density. Alternatively, 
one could argue in the following way: the phase fluctuation term 5A k on the r.h.s. of (20) 
can be thought of having its origin in a replacement of the scalar (<3> ext ) and vector (A cxt ) 
potentials in the first term on the r.h.s. of Eq. (20), representing the external potential 
contributions to 5e k and 5e k , by their gauge-invariant counterparts 

/F.cxt ^ /IS ext 

ldx 
cdt' 

A cxt _^ A cxt + _ ( 22 ) 

The phase variable \ characterizing this gauge transformation can then be fixed by the 

requirement of charge conservation, which, together with the trivial connection between \ 

and the order parameter phase fluctuation 5(f), 

ci 2e 

o<P = X, 

c 

immediately leads to the result (19). This demonstrates clearly the equivalence of the 
existence of the Anderson-Bogoliubov mode with gauge invariance as well as the connection 
of gauge-invariance with the conservation law for the (charge) density. 

Next we would like to demonstrate two simple consequences of a gauge-invariant for- 
mulation of superconducting response theory, namely those for the current response in the 
static {uj — > 0) and for the Raman response in the homogeneous (q — > 0) limit. The static 
limit of the kinetic equation (20) reads 

dn k - —— be k = \ k be k - \ k r) k — —. (23) 
oik E P \% 

Integrating this according to the prescription (16) one gets the static (super) current response 

expressed in a generalized London formula: 

je(q, 0) = ej(q, 0) = -— A (q, 0)A cxt (q, 0), 

c 

Aiq.^-j A A q : q - A ) : (21) 

{ q- A q 

A (q, u) = 2 A P v p : v p , 
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Clearly, the second term in (24), sometimes referred to as backflow term ||20|| , is neces- 



sary to maintain charge conservation in the general case of an anisotropic superconductor. 
Equivalently, it guarantees that the current is purely transverse in the static limit. 

Let us now turn to the Raman response in the homogeneous limit. Here we may ignore 
terms linear in <3> ext and A ext . For q — * the kinetic equation (20) assumes then the 
strikingly simple form 

5n k = -X k 5ei + (25) 

in which the second term on the r.h.s. originates from the order parameter phase fluctuations. 
The Coulomb interaction becomes irrelevant in this limit, as will become clear later, and we 
may write Se k = r y k 5e 1 . The homogeneous Raman response of anisotropic superconductors 
can then be written in a form analogous to the London limit of the current response as 

<5ra 7 (0,co>) = A 77 (0,u;)<5e 7 (0,c<j), 

A afe (q, u) = - ^X ab - Xa ^ b | (q, w), (26) 

A afe (q,^) = 2^ A p a p b p ; a p ,b p = l,q p . 
p 

As in the case of the current response there is a "backflow" term (the second in the curly 
brackets), which guarantees charge conservation. This is easily seen in the limit of parabolic 
bands, where the Raman vertex 7^ is a constant and as a consequence the two terms in 
curly brackets cancel precisely. This is just another way of stating that there are no density 
fluctuations possible in the homogeneous limit q — > of a superconductor. 

It is worth noting that there exists a homogeneous limit of the electronic Raman response 
also in normal metals when quasiparticle scattering processes are important characterized by 
a momentum-dependent lifetime 17. In Appendix B we show that an equation similar to (26) 
holds in the normal state with the Tsuneto function X k replaced by A k = (— dn° k / d^ k ){l — 

Let us now turn to a description of the Raman response at finite wavevector q. Our 
starting point will be Eq. (8) in which we select the even-parity contributions 5n^ and 
5t k = 5e + 7fc5e 7 with 5e = V q 8ni\ 



5nt = [S k + n 2 X k }5e + [ lk S k + fi 2 A fc ^]<5e 



Ai 



u 7 , 



bk ~ 2 2 — ' 
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Integration over momenta k in Eq. (27) yields equations that describe the coupling of density 
and Raman response: 

6n, = x^Se + x^Se,, (28) 
where we have defined generalized (free) response functions 



co — u q An p [ a; - ?? p 



The quantities (l; u ) are straightforward generalizations of the free superconducting 
Lindhard function which include vertices a^, bf.. The Anderson-Bogoliubov collective mode 
causes the second term in the transverse contribution A a b (note that A a i =Aj&= An = 0) 
and the longitudinal term which is characterized by the gauge mode frequency u q . Finally 
we explicitely work out the Coulomb renormalization 5eo = V q 5rii which couples Raman and 
density fluctuations and arrive at the following final result: 



7' 



, ^ A/ (0 M 0) V (0) V (0) 1 

. .(0) Xal Xlb , Xal Xlb 1 /o n s 

Xab ~ Xab (0) 1 (Of- -> \ 6() ) 

Xu Xu 
e = l-V g x?l 

It should be noted that the Raman response in superconductors in the small q limit in a form 
equivalent to Eq. (30) has first been derived by Klein and Dierker [^TJ using the Greens- 
function method. Our result from the kinetic equation method looks slightly different and, 
particularly, one can see that the Coulomb interaction acts so as to split the Raman response 
into an unscreened "transverse" and a dielectrically screened "longitudinal" part, the latter 
being described by the full dielectric function of e(q, u) of the superconductor in complete 
analogy to the behavior of the current response, discussed in Ref. [p0 |. Furthermore, an 



inspection of Eq. (30) shows, that all terms in the full response function Xab except A 77 
are at least of the order 0(q 2 ), the last (longitudinal) term being even smaller, of the 
order 0(q 2 /e). The role of the Coulomb interaction is thus limited to show up in terms 
of the order of at most 0(q 2 ) and is seen to be negligible in the homogeneous limit. The 
role of the collective Anderson-Bogoliubov mode, on the other hand, mainly consists of 
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providing a particle number conservation law. This manifests itself first in the correct mass 
fluctuation limit of Eq. (30), in which 7fc=const and all transverse terms vanish identically 
leaving essentially the screened RPA Lindhard response of the superconductor in which the 
collective mode gets shifted to the plasma frequency. Second it manifests itself in "partial 
screening" effects in the homogeneous limit of the Raman response, described by the second 
term of A~~. These effects will be discussed in detail later. 



Since g£ << 1 in the cuprates (with £ the coherence length), we are mostly interested 
in Raman scattering with vanishing momentum transfers. For such a case it is essential 
to conclude at this stage that the most important contribution to the electronic Raman 
effect in superconductors in the small-q collisionless limit comes from the response function 
A 77 (0, uj) = lim q ^ X77(q? w). Physically, this function describes the photon-induced break- 
ing of a Cooper pair into a pair of Bogoliubov quasiparticles with total momentum q which 
is approximately zero. 

Writing the Raman vertex 7 as a sum 7(k) = 7 +<57(k) of an isotropic and an anisotropic 
part (using Fermi surface harmonics 5j(k) = J2l^o1l(X)), one may decompose the screened 



Raman response function at zero temperature in the limit of small momentum transfers as 
X77 = X|| + X± with 



into a longitudinal part (||), affected by (longitudinal) screening through the dielectric func- 
tion e of the superconductor, and a transverse part (_L), independent of e. Thus for q — > 0, 
only the transverse part of x remains. We also see that in the case of the isotropic density 
vertex, 7(k) = 70, the long range Coulomb forces completely screen the response and thus the 
only contribution to Raman scattering at q = comes from energy bands with nonparabolic 
dispersion, i.e., the L = 2 and higher terms, representing intracell charge fluctuations. 

Taking the limit of q — > and carrying out the £ integration in the Tsuneto function, 
we obtain the final result for the Raman response at finite temperatures, 



B. Final Results for q — > 0. 




(31) 



(32) 
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where the spectrum of Xa,b is given by 



*-W & /° (k)i(k)|AM 



V4T 



4 I Afk) I 2 , 



(33) 



where Np is the density of states for both spin projections, Re denotes taking the real part, 
and (• • •) denotes performing an average over the Fermi surface defined by 

Jd 2 k5{E F -e{k))A{k) 



J d 2 k5{E F - e(k)) 



The real part of x is given as 



V:,,^') = (a(k)/;(k).4'(^))tanl,(^) 



(34) 



A'(u) = N f 



A(k) | 2 J 



arctan 



v/(^/2)2-|A(k)|- 



log 



2 v /|A(k)P-(^/2)2 
W2-y/(W2) 2 -|A(k)P 



W /2 +v /( W /2) 2 -|A(k)|2 



A(k) | 2 > (^/2) 2 , 

A(k) | 2 < (cu/2) 2 . 



This is the expression for the gauge invariant Raman response which has Coulomb screening 
and the Anderson-Bogoliubov gauge mode taken into account. It is also the form for the 
response calculated diagrammatically in the "pair approximation" for the bare bubble, mod- 
ified with the usual R.P.A. screening terms |^T],^2j. It does not take into account any vertex 
corrections resulting from the pairing interaction in other channels other than the pairing 
channel, as explained in the preceeding discussion. We will refer back to this expression in 
the following Sections. 

Important information can also be obtained from the temperature dependence of the 
response in the limit of zero frequency shifts, i.e., the static reponse. It can be shown that 
the ratio of the response in a superconductor to that of a normal metal in the limit of 
vanishing frequencies is given by the simple expression 



0) 2(/(| A(k) |) | 5 7 (k) | 2 ) 



(35) 



X".s.(^0) (|57(k)l 2 ) 

where / is a Fermi function p3| . This is an exact result which does not depend on vertex 
corrections and only depends appreciably on impurity scattering for nearly resonant impurity 



scatterers 23,24 



The important feature in all these expressions is that in general a coupling of the Raman 
vertex to the energy gap can occur under the momentum averaging over the Fermi surface. In 
all previous studies, this fc-dependence was either ignored or not fully exploited to determine 
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important information concerning the symmetry of the energy gap. This will be explicitly 
demonstrated in the following section where we evaluate these expressions for the screened 
Raman response for various pair state candidates and discuss its relevance to the cuprate 
materials. 

We close this section by returning to the question of the presence of massive modes and 
final-state vertex corrections (see above). In principle, the massive modes can lie at low 
frequencies and affect the low frequency behavior of correlation functions, and in particular, 
could even be used as a signal for a certain type of order parameter symmetry. We have 
carried out an analysis of the full gauge invariant calculation in the Appendix C for both 
cylindrical and spherical Fermi surfaces for a generalized pairing interaction. We identified 
the position of the collective modes as a function of coupling strength for a gap of d x 2_ y 2 (T3 
representation using the notation of Sigrist and Rice, |[25|| ) symmetry for both a spherical 
and cylindrical Fermi surface. Our conclusions are threefold: 1) the Goldstone modes do 
affect the Raman spectrum in the limit of q — > 0, in certain polarization symmetries (Ai g ) 
where the "backflow term" in (32) is finite, 2) optical (massive) collective modes couple to a 
light probe only for case of a Fermi surface with ^-dispersion for <i-wave vertex corrections 
in the B\ g and E g channels, 3) these modes are significantly damped and have a vanishing 
residue for those modes which lie at low energies. Therefore, the collective modes are of 
little importance to the Raman spectrum. These conclusions can be generalized to other 
energy gaps which have line nodes on the Fermi surface. Also, we discuss the role of the 
final-state interactions and find that while the overall shape of the spectra can be affected, 
the corrections are relatively minor. The details are contained in Appendix C, including a 
more general discussion of the role of vertex corrections. Therefore, we can conclude that 
Eqs. (32-35) give an adequate desciption of the Raman response. 



III. BAND STRUCTURE, FERMI SURFACE AND THE RAMAN VERTEX FOR 

TETRAGONAL SYMMETRY 

In this section we aim at providing a link between the Raman vertex and band structure 
for tetragonal crystals. In particular, we show how the choice of light polarization orienta- 
tions results in selection rules for the symmetry components of the Raman tensor. In the 
first subsection we only consider scattering within a single band while we consider the case 
of multiple bands at the Fermi surface in the following subsection. 
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A. Single Band 



As we have seen in the previous section (see Eq. (2)), the Raman tensor is directly 
related to the curvature of the band dispersion. In the following we limit our considerations 
to tetragonal materials which are relevant to high T c superconductors. Although these 
materials possess orthorhombic distortions away from tetragonality, the selection rules for 
example from phonon scattering in the normal state seem to indicate that these orthorhombic 
distortions are small and that a tetragonal symmetry classification can be employed with 
little loss of generality. A simplest choice for the band structure which contains the basic 
physics of 2D-like tetragonal systems with lattice constant a is given by, 

e(k) = —It [cos^a) + cos(k y a)] + At 1 cos(k x a) cos(k y a). (36) 

Here t and t' are the nearest and next nearest neighbor hopping parameters, respectively. 



This is the anti-bonding band derived from a reduction of a three band model |26[ , which 
gives by far the largest contribution to the density of states at the Fermi level for the cuprate 
systems |[27|| . The 2D Fermi surface is defined through the relation e(k) = fi, where /i is the 
chemical potential, which in turn determines the Fermi momentum, 

V sm (w / 

The scalar prefactor kp{<p) can be expanded with respect to the fully symmetric basis func- 
tions (A\g or using the notation of Sigrist and Rice [25]) for the tetragonal D 4h point 
group, 

k F ((p) = kf + kf cos(4<^) + 4 4) cos(8<^) + ■ • • 

= kf+ £ k^aKv). (38) 

L=2,even 

The higher order Fourier coefficients k^ for L > take into account deviations from 
cylindricity of the Fermi surface, kp\ The basis functions for the irreducible representations 
of the point group symmetry can be generalized as 

ai(p) = cos(22^+^), JG{0,1}, (39) 

where a J L transforms according to the A,- +l9 -symmetry, and 

jn 



b 3 L (ip) = cos 



(2£ - 2)<p , 2 
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je{0,i}, (40) 



which transform according to the 5j + i 9 -symmetry . 
The Raman tensor is given by 

7(k)= mAf _1 (k) 

^ 2 ( cos(kF X a)[t — 2t' cos(kF y a)] 2t' sin(kp x a) sm(kF y a) \ 
y 2t' sin(/ci^a) sin(/cp y a) cos(fci^a)[t — 2t' cos(/cF X a)] / 

where = de(k)/dk a dkp. We denote matrices in fc— space (except the Pauli matrices) 
with a double-arrow. We proceed to expand the Raman tensor in quarternions, 

3 



-l A 



7(k) = mM (k) = ^ (m) (k) r r 

m=0 



(42) 



with r m the Pauli matrices in 2-D k-space. The polarization light vectors select elements 
of the Raman tensor according to 



7 (k) = e'7(k)e 5 = ^(eVe 5 ) <5 m (k). 

m 

Next we expand the functions 5^ (k) in terms of Fermi surface harmonics 

L=2,even 

where 



(43) 



(44) 



(45) 



'a£foO> A la (m = 0), 

B 2s (m = l), 
(m = 2), 
[b°M, B lg (m = 3). 

The coefficients eV m e of the Raman vertex can only assume the values or 1 and have the 
following explicit form: 



e l r m e s 





+ ejej], 


A 9 


(m 


= 0), 








(m 


= 1), 


Fx ?/ 


- eVl 




(m 


= 2), 


, [ e x e x 


- eVl 


Big 


(m 


= 3). 



(46) 



We now can reconstruct the Raman vertex (see Eq. (2)) as 

7(k) = E7rCM 

L,m 

7 ™ = eV m e s 4 m) . 



(47) 
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We therefore see the connection between polarization orientation and symmetry. For ex- 
ample, if e 1 = e s = x are chosen, then a combination of A\ a and B\ g is selected. If 
e 1 = e s = 2~ 1 / 2 (x + y) then A\ g and Bi g symmetries are combined. These examples show 
in particular that A\ g symmetry cannot be individually accessed for in-plane polarizations, 
and subtraction procedures must be used. 



B. Multi sheeted Fermi surfaces 



The above consideration can be easily adopted to the case of multiple energy bands 
contributing to the Fermi surface, for instance, due to the contribution of the chains as 



well as the planes in YBCO |27]. The total electronic Raman cross section for intraband 
scattering is just the sum of the contribution to the cross section for each band and thus 
results from the addition of the Raman response functions, 

Xtotal = Y.Xa{&a}, (48) 

a 

where Xa{A a } denotes the contribution to the scattering from band a with an energy gap 
A a for each band. The Raman vertex for each band is separately calculable in the same 
manner as in the previous section, and the same considerations there can be carried over to 
the multi-sheeted case trivially. 

We draw the important point however that the contribution from each band will be 
weighted by the relative density of states and curvature of each band at the Fermi level. 
Therefore, for the case of the cuprates we believe that by far the largest contribution to 
electronic Raman scattering is thus arising from the single anti-bonding Cu-0 layer band 
with the largest density of states and greatest curvature at the Fermi level. Further, it is 
not even clear whether an energy gap exists on the other energy bands (e.g., associated 
with the chains in YBCO), and if it does exist what symmetry it has. Therefore, we feel 
that to a very good approximation, the intraband scattering has a dominant contribution 
from the anti-bonding band. This is certainly the case for the single Cu-0 layer compound, 
Tl 2 Ba 2 Cu0 6 . 

We close this subsection with a discussion concerning interband scattering, or scattering 
from the anti-bonding band to another band near or at the Fermi level. In this case a 
separate theory is necessary to take into account the separate band dispersion and energy 
difference of the two bands, and the energy gap on each of the bands. This at present remains 
a future consideration. However, we remark at this point that once again this contribution 
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to the scattering cross section will be smaller than the contribution from scattering within 
the anti-bonding band due to the small curvature of the other energy bands and the value of 
the density of states at the Fermi level (again, the chains for example). Also, since the nature 
of the energy gap on the chains is completely unclear, we consider interband scattering to 
be beyond the scope of the present manuscript [ f28|j . 



IV. PREDICTIONS FOR VARIOUS PAIR STATES 



In this Section we evaluate Eqs. (32) and (33) for various pair states that are discussed 
in the literature as candidates to describe the pairing symmmetry in the cuprates. We will 
first discuss the simple case of an isotropic gap and the discuss ci-wave gaps (in particular, 
we will focus on the d x 2_ y 2 pairing symmetry) and then discuss gaps which are anisotropic 
but do not contain line nodes, e.g., an s + id and an anisotropic s-wave gap. 



A. Isotropic s wave 



For the case of an angular independent gap, A(k) = A , the Raman response is given 
by the simple BCS expression. It has the stereotypic features of BCS theory, namely, 
the existence of an energy gap threshold 2A required to break a Cooper pair, and the 
ubiquitous square root divergence associated with the gap edge. We note in particular 
that the Raman vertex couples trivially to the energy gap and just determines an overall 
prefactor governing the Raman intensity. The vertex does not affect the lineshape and thus 
the spectrum is polarization independent apart from an overall prefactor. A polarization 
dependence can be generated in BCS theory by taking into account channel-dependent final- 
state interactions P]] and/or impurity scattering [2~3~|, and accurate fits to the Raman data 
on A- 15 superconductors can be obtained. However, for the most part this only produces 
an channel dependence in the vicinity of the gap edge and thus the main feature of the 
response is the uniform gap existing for all polarizations, which clearly cannot give an 
adequate description to the Raman spectra of the cuprate materials. 



B. d x 2_ y 2 pairing 

The explicit symmetry dependence of the spectra results due to a coupling of the energy 
gap and vertex when the energy gap is anisotropic. To demonstate this, we now carry out 
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the averages over a cylindrical Fermi surface in Eq. (33) using a d x 2_ y 2 gap, A(k, T) = 
A (T) cos(2<y?) (see Appendix A). We find that the Raman spectrum can be written down 
analytically in terms of complete elliptical integrals. Taking screening into account and 
defining x = u/2A , we obtain for T = 0, 

x'iC =X , B 1 ,(q = 0,u;)= (49) 
2^F7l lg x f [(2 + x 2 )K(x) - 2(1 + x 2 )£(x)], x < 1 , 

3vra; X + 2x 2 )fT(l/x) - 2(1 + x 2 )£(l/x)], x > 1, 

xb:: = xb 2 m=°^)= ( 5 °) 

2jV F7 | 29 x f [(1 - * 2 )K(x) - (1 - 2x 2 )E(x)}, x<l, 
3ttx X \ ;r[(2 - 2x 2 )K(1/x) - (1 - 2x 2 )E(l/x)}, x > 1, 

i.e., the _Bi 9 and i?2 9 channels are not affected by Coulomb screening. This is consistent with 
the mass fluctuations being only intracell in nature for these symmetry channels. However, 
the Ai g channel which contains both inter-and intracell fluctuations is partially screened 
and is determined via 

xt^)=XA lg ^)-^^, (51) 
with the spectral functions 

(q = 0,o;) = ^^ (52) 
!<" is 157rx 



X 



[(7 - 8x 2 + 16x 4 )K(x) - (7 - 12a; 2 + 32x 4 )£(x)], x < 1, 

x[(U - 28x 2 + 32x 4 )K{l/x) - (7 - 12x 2 + 32x 4 )E{l/x)}, x > 1, 



„ . n , v^iV F7Al r[(l + 2a; 2 )^(a;)-(l + 4 a ; 2 ) J E(a;)], x < 1, 

^4, 1 q = 0,w = < 53 

191 3vrx \x[(-1+x 2 )K(1/x)-(1 + x 2 )£(1/x)], x > 1, 

and 

y, , q = 0,w = < 54 

A1 ' 1VH ^ 7TX \x[X(l/x) -^(l/x)], X>1. ^ ' 

The real parts are obtainable through Kramers-Kronig transformation or numerically inte- 
grating Eq. (34). The response functions for finite T are obtained simply by multiplying Eqs. 
(49-54) by the factor tanh(cj/4T). The partial screening of the A lg channel arises technically 
from the observation that the square of the energy gap enters into the response function in 
Eq. (33). For the case of d x i_ y i pairing symmetry, the energy gap squared contains a term 
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which transforms according to A\ g symmetry which leads to a finite overlap with the A\ g 
vertex in Eq. (32). This corresponds to partial "transverse screening" of the A\ g channel, 
and via this mechanism the intercell fluctuations are removed. Similar considerations hold 
if the gap were of another <i-wave symmetry other than d x 2_ y 2. 

These functions are plotted in Fig. (1) for the 3 symmetries indicated. We immediately 
see that the spectra is extremely polarization dependent, in contrast to the case of isotropic 
s-wave superconductors which is dominated by the square root divergence at the threshold 
in each channel. We see that the peak in the Raman spectra lies at different frequencies 
Upeak ~ 2A (T), 1.6A (T) and 1.2A (T) for the Bi g , B 2g and Ai g channels, respectively [|30| . 



The symmetry dependence of the spectra is a direct consequence of the angular averaging 
which couples the gap and Raman vertex, and leads to constructive (destructive) interference 
under averaging if the vertex and the gap have the same (different) symmetry. Thus it has 
been reasoned that the symmetry which shows the highest peak position gives an unique 
indication of the predominant symmetry of the gap 0. The peak positions can be mildly 
affected by including interaction vertex corrections as discussed in Appendix C, with the 
net result being a slight upward shift of the peak location in the B 2g and A\ g channels. 
Also, while the presence of ^-dispersion has little effect on the B lg (apart from cutting off 
the logarithmic divergence) and B 2g channels, the A lg peak position can be changed due to 
the addition of a terms which has its main contribution at slightly lower frequencies. This 
effect is small provided that the Fermi surface is mostly cylindrical. We refer the reader to 
Appendix C for details. 

The symmetry dependence is also manifest in the low frequency behavior, which can be 
written as 

XB ls (^-0) = 3iV F 7l l9 X 3 /4 + O(x 5 ), 
XB 2g (u^0) = N Fl 2 B2g x/2 + O(x 3 ), 

A> - °) = N Fl 2 M x/2 + 0(x 3 ), (55) 

i.e., the spectrum rises slower in the Bi g channel than the A lg or B 2g channels, which have 
the same linear rise with frequency. The power-laws are insensitive to vertex corrections 
and arise solely due to topology arguments. The appearance of power laws are a signature 
of an energy gap which vanishes on lines (points in 2D) on the Fermi surface. However, 
the channel dependence of the exponents are unique to a d x 2_ y 2 pair state. These channel- 
dependent power laws have been observed in the electronic contribution to Raman scattering 
in BSCCO |f7|,|12|, YBCO [Tl|] and TBCO |13| , |3T| and are strong evidence for a d-wave gap 
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of this symmetry as opposed to d xy , d xz or d yz symmetry, which also have nodes on lines on 
the Fermi surface |]7|j30|. 

We now investigate the temperature dependence of the theory and contrast to that of 
an isotropic s-wave superconductor. Using a weak coupling expression for the temperature 
dependence of the energy gap (2A /T C = 4.2794) (see Appendix A), we numerically evalu- 
ate Eq. (35) for the temperature dependence of the normalized static response while taking 
screening into account. This function describes how the gap in the Raman spectrum at low 
energies opens up with cooling below T c . The results are plotted in Fig. (2) as a function 
of T/T c for a d x 2_ y 2 energy gap compared to a BCS isotropic gap. The low temperature be- 
havior is given by a power-law in T for all channels for the d-wave case while the ubiquitous 
exponential dependence in T is seen for all channels in the s-wave case. The power-law 
behavior for the ci-wave case is channel dependent, with exponents identical to those of Eq. 
(55), in the sense that u can be replaced by T. What is remarkable is that the fall off of 
the Fermi function at low temperatures is quite slow in those channels which are orthogonal 
to the symmetry of the gap, where the A\ a and B2 g channels show a residual broadening 
at T/T c = 0.3 of roughly 20 percent of that of the normal state. This was argued in the 
case of electronic Raman scattering to be further evidence for an energy gap in the cuprate 
materials which has predominantly B lg character, due to the observation that a gap opens 
up quickly in the Bi g channel compared to others channels which have been probed via 
Raman scattering 0. 

As we have remarked before, since the Raman density response function only depends on 
the magnitude of the energy gap, it cannot be sensitive to the phase of the order parameter 
and thus cannot be directly used to determine if the gap changes sign around the Fermi 
surface. While the power-law behavior at low frequencies and/or temperatures is indicative 
to the presence of nodes, in principle a very highly anisotropic order parameter with a 
small uniform gap everywhere on the Fermi surface could mimic the behavior of a gap with 
nodes when inelastic scattering or experimental resolution smears the threshold. In principle 
the detection of a threshold can then only be performed at very low temperatures where 
activated behavior can be observed. Since this remains a possibility we now discuss two 
types of energy gaps which are anisotropic but have a finite gap around the Fermi surface. 
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C. s + id pairing 



Kotliar and Joynt have suggested the possibility that an order parameter which is a 
superposition of an s-wave and a d-wave gap can also provide an adequate description to 
the various transport and thermodynamic measurements on the cuprate systems f32[. The 
s + id state has the interesting feature that in pure tetragonal superconductors there would 
be two transition temperatures associated with the formation of each gap separately. This is 
a consequence of the admixture of two different representations of the energy gap. However, 
orthorhombic distortions remove the x <-> y symmetry and thus A\ g and B\ g belong to the 
same representation. This leads then to one transition temperature. Since the orthorhombic 
distortions are quite small in the cuprates (judging from the observed phonons and selection 
rules), the transition temperature will be broad, which is in some conflict with the resistive 
transitions seen in the cuprates. 

Nevertheless, we investigate what such a gap predicts for the Raman response by 
evaluating Eq. (33) on a cylindrical Fermi surface using a gap of the form A(k) = 
A S {T) + zA^(T) cos(2<£>). The results can written again in an analytic form in terms of com- 
plete elliptical integrals. Taking screening into account and defining x 2 = [(u/2) 2 — A^]/A^, 
we obtain for T = 0, 




(56) 



x 



( 



[(2 + x 2 + 3A 2 /A 2 )K(x) - (2 + 2x 2 + 3A 2 J A 2 d )E(x% x<l, 
x[(l + 2x 2 + 3A s /A 2 )K(l/x) - (2 + 2x 2 + 3A 2 J A 2 )E(l/x)}, x > 1, 




(57) 




(58) 



with the spectral functions 




(59) 
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f [(7 - 8x 2 + 16x 4 - 5A 2 /A 2 (1 - Ax 2 ))K(x)- 
(7 - 12x 2 + 32x 4 - 20A 2 /A 2 (1 - 2x 2 ))E(x)], 
' ^ x[(23 - 40x 2 + 32x 4 + J£%(3 - 8x 2 + 8x 4 ))K(l/x)- 



(11 - 28x 2 + 32x 4 - §&(1 - 2x 2 ))£(l/x)], 



x < 1, 



X > 1, 



and 



Xa 



ly 



,i(q = o,w) 



[(i 



X < 



x "(q = 0,^) = e(x 2 ) 



37ra; 

2x 2 + 3A 2 /A 2 )K(x)- 
{1 + Ax 2 + QA 2 /A 2 d )E(x)], x<l, 
-x[{\ - Ax 2 + - 2x 2 ))K(l/x) + 

(l + 4x 2 + 6A 2 /A 2 )£(l/x)], x>l, 



2N F A d { [(1 + A 2 /A 2 )K(x) - £(x)], x < 1 , 
ttlu \x[{l + ^)K{l/x)-E{l/x)}, X>1. 



(60) 



(61) 



The results are plotted in Fig. (3) using an value of A s /A d = 0.25. The finite gap A s is 
responsible for the threshold appearing at 2A S , which is the minimum energy required to 
break a Cooper pair. The spectra are polarization dependent, with the B 2g and A lg spectra 
displaying a discontinuous jump at the threshold while the B\ g channel shows a continuous 

= 2wA 2 + A 2 ,. Since both the A lg and B 2g channels 



rise from zero to a peak at to = 2A ma 
exhibit broad maxima for the case of a pure d x 2_ y 2 gap, the presence of the 2A S threshold 
will imply a shifting of the peak of the spectra away from the pure case to values lower in 
frequency if 2A S < u pea k (or at least create a shoulder at 2A S ), while also removing any 
difference in the peak position between the A\ g and B 2g channels. Taken by itself this cannot 
be reconciled with the experimental data unless of course A s is very small. 



D. Anisotropic s wave pairing 

Starting from band structure arguments, recently an anisotropic s-wave energy gap of 
the form 



A(k) = A + Ai cos 4 (2(^) 



(62) 



has been proposed also to explain the cuprate materials |33]| . This energy gap is anisotropic 
with predominantly I? lg -like character but does not change sign around the Fermi surface. 
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This is one specific example of an anisotropic s— wave energy gap. While other represen- 
tations for anisotropic s— wave gaps of course do exist (in particular, a possibility recently 
suggested by photoemission is A(k) = cos(k x a) + cos(k y a) = const. + Acos(4y9) oc cos 2 (2</>) 



34j), we remark that the response calculated here is not qualitatively different from other 



cases and thus address only this one case. 

Again we evaluate the Raman response for such a superconductor numerically and plot 
our results for Ao/Ai = 0.25 in Fig. (4). Immediately we see similar behavior as in the 
previous case with one notable exception. While the spectra each show a 2Ao threshold, the 
B2 9 channel displays a B.C.S.-like singularity at the threshold and the A\ g a large increase 
near the threshold that removes any trace of a peak-like structure in the spectra at higher 
frequencies. Again this would predict the same peak position (or a shoulder) for the A\ g and 
B2 9 channels. Therefore, for Ao not too small, it is not possible using just the symmetry of 
the gap alone to arrive at a situation where the peaks in the Raman spectra in the A\ g and 
B2 9 channels lie at separate high energies, again which is not in agreement with experiments. 

In principle the s-wave component of the gap can be made vanishingly small. However, 
the greater anisotropy of the energy gap compared to the d x 2_ y 2 case leads to a further 
anisotropy of the position of the peaks in each channel, with peaks in the spectra at uo = 
2A max , 0.6A max , and 0.2A max for the B lg , B 2g , and A lg channels, respectively, where A max = 
A + Ai. On top of this, the low frequency power-law behavior is linear in each channel, in 
contrast to the d x 2_ y 2 results. This is a general feature that the spectra become more and 
more polarization dependent the greater the anisotropy of the energy gap (compare to Fig. 
(1) for the d x 2_ y 2 case). 



V. COMPARISON WITH DATA ON THE CUPRATE SYSTEMS AND 

CONCLUSIONS 

In this Section we present a comparison of the theory for a d x 2_ y 2 paired superconductor 
to recent measurements on the electronic Raman continuum in three cuprate superconduc- 
tors. In what follows, we plot S(u), Eq. (1), which is given by x of Eqs. (49)- (54) mutliplied 
by the Bose factor. In drawing the fits to the spectra, the following procedure is employed. 
First the fit to the B\ g spectrum is made which determines the maximum value of the 
energy gap Ao via the position of the peak in the spectrum. Next, the derived response 
is convoluted with a Gaussian which mimics the effect of finite z-dispersion of the Fermi 
surface, experimental resolution, inelastic scattering, etc. Once this is done, the parameters 
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remain fixed and only the prefactor of the vertex is left to be adjusted to match the overall 
intensity, which has no effect on the lineshape. Since 7 is in principle derivable from band 
structure but presently unknown even for such simple metals as Aluminum, this remains a 
free parameter. 

We first fit the data taken on single crystals of as grown Bi 2 Sr 2 CaCu208 (T c = 90K) 



obtained in Reference |L2[ for all symmetries at T = 20K, where a subtraction procedure 
has been employed to ascertain the A\ g signal (see Section 2.2). The comparison of the 
theory with experiment is shown in Fig. (5). The parameters used to obtain the best fit to 
the spectrum are Ao = 287 cm -1 and a smearing width of T/Ao = 0.15. The theory gives 
good agreement with the data at low frequency shifts while at higher frequency shifts the 
theory fails to produce the broad continuum which is relatively constant up to the scale of an 
eV. This is most likely due to the neglect of impurities and/or electron-electron scattering 



135|, which is beyond the scope of the paper. We see immediately that the peak positions 
in the B 2g and A\ g channels given by the theory automatially agree with the data. Also 
the asymptotic behavior of the continuum at low frequencies given by Eq. (33) is shown 
in the data when one neglects the phonons at roughly 100 and 330 wavenumbers. Again, 
these power-laws are intrinsic to a d x 2_ y 2 pair state. Lastly, the ratio of the intensity of 
the spectra in different channels is consistent with Eq. (35) and Fig. (2), which predicts 
that the Bi g channels shows the smallest intensity at low frequencies while the Ai g channels 
shows the largest. All of the experimental features are thus consistent with the theory at 
least at low frequency shifts u < 1000 cm -1 . 

We now turn to the data taken on single crystals of YBa 2 Cu30y (T c = 88K) obtained 
in Ref. [J l.| for all symmetries at T = 20i^, where the same subtraction procedure used in 
BSCCO was employed to ascertain the A\ g signal. The comparison of the theory to the data 
is shown in Fig. (6), with the parameters A = 210 cm -1 and r/A = 0.2 Again, the theory 
gives a good description of the data for uj < 1000 cm^ 1 . We again see the peak positions at 
relatively the same place as in BSCCO, and power-laws linear in frequency at low shifts for 
the B 2g and Ai g channels. We note that while the cubic rise of the spectra predicted by the 
theory fits rather well with the Bi g data in BSCCO, we remark that the Fano effect of the 
Bi g phonon which appears to be stronger in YBCO than in BSCCO can obscure the rise of 
the spectra at low frequencies at give the appearance of a linear dependence on frequency 
|3l]fl . Lastly, the ratios of the response in the static limit again are consistent with the theory. 

Lastly, we investigate the single layer Thallium compound Tl 2 Ba 2 CuOg (T c = 80K) 
obtained in Ref. [|13||. The sample is most likely the most affected by disorder and therefore 
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our theory will not be expected to give the best fit to the data. Our fits is shown in Fig. 
7 for the B\ g and the mixed A\ g + B2 g channels at T = 20K, where we have ignored 
the contribution from the B2 g channel since it is believed to be minor compared to A\ g . 
The value of the parameters used are Ao = 232cm _1 and smearing width T/Aq = 0.25. 
All phonons have been subtracted. Once again the theory gives a good description of the 
relative peak positions. Considering also that this compound has only one Cu-0 layer, the 
agreement of the theory with experiment also validates the assumption that the Raman 
scattering results predominantly from intraband fluctuations of the single Cu-0 layer band, 
and that interband scattering can be neglected. The theory also accounts for the linear rise 
of the spectrum for the A\ g channel for low frequencies. However, the theory cannot account 
for the linear rise of the spectrum in the B\ g channel. This most likely is due to the neglect 
of impurity scattering. The residual scattering near zero frequency shifts is also borne out 
by the theory, although the amount is underestimated for the B\ g channel. Again, this most 
likely has to do with the neglect of impurity scattering. 

We remark that the theory is incomplete in that the theory fails to describe the flat 
continuum at large frequency shifts. Moreover, the theory cannot be extended to the normal 
state since the response functions vanish at T c in the limit q — > due to phase space 
restrictions. Here the additional physics of electron-electron scattering and/or impurity 
scattering must be incorporated to have a consistent theory to simultaneously describe 
the normal and superconducting state data. This remains a topic of further research |24| . 
Nevertheless, the low frequency behavior of the spectra, and in particular, the relative peak 
positions of each polarization channel are quantitatively described by the theory. 

Thus we have seen that the Raman measurements on the cuprate systems provide a large 
body of symmetry dependent information all of which agrees with the predictions of d x 2_ y 2 
pairing. Of course at present the information from Raman alone cannot completely rule out 
the possibility of the presence of a very small gap which exists over the entire Fermi surface 
nor can it determine whether the gap changes sign around the Fermi surface. However, the 
theoretical comparison shows that the gap must be predominantly of B\ g character, and the 
low temperature and low frequency data seem to indicate that the s-wave component of the 
gap must be very small if it exists at all. More precise measurements could of course clarify 
this point further. Also, more work is needed from band structure to pin down magnitude 
of the Raman tensor elements to predict the overall intensities. 



26 



ACKNOWLEDGMENTS 



The authors would like to thank Drs. R. Hackl, G. Krug, R. Nemetschek and B. Stadlober 
for providing us with their data and discussions. Similarly, we thank D. Reznik, R. T. 
Scalettar, A. Virosztek, A. Zawadowski, and G. T. Zimanyi for enlightening discussions. 
This work was supported in part by N.S.F. Grant number 92-06023 and by the American 
Hungarian Joint Grant Number NSF 265/92b. One of the authors (T.P.D.) would like to 
acknowledge the hospitality of the Walther Meissner Institute, the Research Institute for 
Solid State Physics of the Hungarian Academy of Science, and the Institute of Physics of 
the Technical University of Budapest where parts of this work were completed. 



APPENDIX A: WEAK COUPLING RESULTS 



In this appendix we want to show that a gap anisotropy causes both the gap at zero 
temperature A (0) and the specific heat discontinuity at the transition AC/Cn (related 
to the slope of the gap function near T c ) to deviate from their respective BCS values of 
(A (0)/T c )bcs = 7rexp(- 7 ) = 1.7638... and (AC/C N ) BC s = 12/7C(3) = 1.4261..., with 
7 = .57721... and £(3) = 1.20205... denoting Euler's constant and Riemann's zeta function 
respectively. A straightforward solution of Eq. (18) leads to 

A (0)_, /A (0)\ 



T, 
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AC 
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2 ln(|A p |/A ) > FS 
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2 >FS 
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N / BCS 
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IA I 2 


^FS 
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IA I 4 


>FS 



(Al) 



These results may be used to generate an interpolation formula for the temperature depen- 
dence of the gap maximum A (T): 



A (T) = 5 sc T c tanh 



7T 



3 AC 



A 2 



S sc \ 2 C N < |A P | 2 > FS VT 



(A2) 



For the special case of a gap with d x 2_ y 2 symmetry, one obtains for the parameters S sc and 
AC/C N to following numbers: 



f expi-7 + if) = 2.5626, 



27rexp (—7 — \ 



12 7 



AC/C.H 7 ^ 5 "°' 6655 ' 



spherical FS, 
2.1397, cylindrical FS, 
spherical FS, 
cylindrical FS. 
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We should emphasize, that these numbers should not be taken too seriously since they 
emerge from a weak coupling treatment. One should rather adopt the convention to treat 
them as parameters which can be adjusted to experiment and so account for strong coupling 
effects in the trivial sense in which they appear as renormalizations of the quantities Ao and 
AC/C N . 

APPENDIX B: IMPURITY SCATTERING IN THE NORMAL STATE 

In this Appendix we would like to demonstrate that the normal state electronic Raman 
response has a nonvanishing q — > limit the presence of quasiparticle collisions. In the 
collisionless limit of a superconductor, the quantity x 77 — A 77 has a finite value in the normal 
state limit A —>■ only at finite wavenumbers q. The normal state Raman response can 
persist in the limit q — > however if one assumes (elastic or inelastic) scattering processes 
to take place. These scattering processes may be generally introduced through a collision 
integral 5I k in the normal state kinetic equation [[1^1 , 



u5n k - q • v fc h k = i5I k {h k }, 
Sh{h k } = —h k + 2 ^ C kp h p , (Bl) 

T k p 

hk = on k - ~^-oe k . 

The structure of the collision integral depends on the scattering mechanism that one is 
interested in and is, in many cases, subject to an approximate treatment of the collision 
operator C kp . However, such approximations have to be consistent with the requirement of 
charge conservation [p7 , p5[ 59[. An approximation with this property is the separate kernel 
approximation fHf . 



The scattering processes are here represented by a momentum-dependent quasiparticle life- 
time r k and a scattering parameter Ao- The continuity equation is obtained from Eq. (Bl) 
by summing over momenta k, 

uSn - q- j = -2iJ2 — (1 - A ). 



p T p 



Particle number conservation is therefore connected with the choice of the scattering pa- 
rameter Ao = 1. One may show that the normal state Raman response in the q — > limit 
assumes form equivalent to the result (26) for the pair breaking Raman effect: 
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A ^0q->0 A77 77 A 



11 



p \ d£ k J 1- tUJT p 

Like in the superconducting case, the projected structure of Eq. (B3) has its origin in 
the particle number conservation law. This result can be applied to impurity or phonon 
scattering as well as to inelastic two-particle scattering and requires the specification of the 
momentum-dependent quasiparticle lifetime for each of these cases. 



APPENDIX C: DIAGRAMMATIC GAUGE INVARIANT RAMAN RESPONSE: 
ROLE OF VERTEX CORRECTIONS AND COLLECTIVE MODES 

In this Appendix we derive expressions for gauge-invariant generalized correlation func- 
tions using a diagrammatic approach with the main emphasis placed on electronic Raman 
scattering. After solving the coupled integral equations for the renormalized vertex, we 
evaluate the position, broadening and residue of the massive and massless collective modes 
as a function of coupling strength for a d x 2_ y 2 energy gap on a spherical Fermi surface. We 
then investigate the collective modes on a cylindrical (2-D) Fermi surface by turning off any 
^-dispersion in the band structure. 

We begin by writing down the expression for the general two-particle response function 
in Nambu space as 

X (q, iw) = ~TJ2J2 T ^[ f (k, q, iu)G(k + |, iu n )i(k, -q)G(k - |, iu n - (CI) 
f(k, q, iu) - 7(k, q) = -T £ ]T r 3 G(p + ^, iu n )t(p, p', iu) (C2) 

iu>n P,p' 

xG(p - iu - iu n )V(p - k + BJlP., k - p + ^~^-)r 3 , 

where Tr denotes taking the trace, Tj are Pauli matrices in Nambu space, and the ver- 
tex 7 determines the correlation function of interest. The dressed vertex T contains the 
interactions V responsible for maintaining gauge invariance. 

Our consideration is focused on the Raman vertex which describes the anisotropic mass 
fluctuations, 7(k, q) = r 3 7(k) with the limit of small q taken. Other choices of the vertices 
are 
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7(k) = r 3 , charge density, 

7(k) = t , spin density, (C3) 
7(k) = kr , current density. 

Thus from Eqs. (C3) we see that while the spin and charge density vertices probe only the 
L = channel and current density probes L = 1, in principle all even L channels contribute 
to the Raman vertex. In fact the charge density is the first term in the expansion of the 
Raman vertex. 

If one replaces the dressed vertex T by the undressed one, Eq. (2), then of course Eq. 
(CI) is manifestly not gauge invariant, and in general the neglect of collective modes arising 
from a gauge invariant treatment could in principle affect the overall spectrum. Usually the 
question of gauge invariance is rather an academic one since the modes that appear in BCS 
systems have little impact on the response functions of a superconductor. It is well known 
that due to the spontaneously broken U(l) gauge symmetry in s-wave superconductors, two 
collective modes appear - an optical one with a frequency of 2A which is damped - and a 
sound-like mode, the Anderson-Bogolubov mode, which is soft and lies in the gap for neutral 
superconductors but is raised to the plasmon energy by the long range Coulomb forces via the 
Higgs mechanism. However, in unconventional superconductors, there can exist in principle 
additional Goldstone modes corresponding to the additional broken continuous symmetries 
such as SO3 spin rotational symmetry in spin-triplet systems plus SO^ orbital rotational 
symmetry in spin-singlet systems if the gap does not possess the full symmetry of the Fermi 
surface. In addition, massive collective modes can arise if the energy gap is degenerate or 
has an admixture of different representations from the point group. The massive modes can 
in principle lie below the gap edge and thus be relevant for the low frequency dynamics of 
correlation functions. For example, the spectrum of collective modes in 3 He is well known 
in both phases and lead to observable effects [TO. However, the collective modes of possible 
d-wave states that might be candidates for strongly correlated systems are not as well 



understood [HJ. There are indications that a ci-wave state of d x 2_ y 2 is particularly favorable 



in systems with strong correlations , thus underscoring the necessity of an understanding 



of the response functions and collective modes for such a superconductor. 

We continue our calculation for the gauge invariant Raman response by first expanding 
the renormalized matrix vertex T along the Fermi surface in terms of crystal harmonics, 

f(k,^) = £f™MC(k), (C4) 
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and do the same for the pairing interaction, 

V(k,p)= £ V^'C(k)^'(p). (C5) 

L,L' ;m,m' 

The integral equation for the renormalized vertex at q = can then be written as 

rp at „ 

frM-7r = — dede' £ (C6) 

Z ia;„ L',L";m',m" 

x (0™;;(k)C;(k)r 3 G(k, za; n - iw)f^(«<;)G(k, zo, n )T 3 ), 

where (• • •) denotes an average over the Fermi surface and N F denotes the density of states 
at the Fermi level. Eq. (C6) is completely general for any type of vertex, interaction and 
gap symmetry. In particular, for the case of an isotropic energy gap, Eq. (C6) recovers the 
previous results for the density, current [J3|, and Raman reponses []21j] . We confine ourselves 



to the case of singlet energy gaps and use the BCS approximation, 

An • \ iu n + £r 3 + A(k)n 

G ^^= (^) 2 -^ 2 (k) ' (C7) 

where E 2 (k) = £ 2 + A 2 (k). More complicated Green's functions could treated with the same 
scheme. However, the analysis gets considerably more complicated and cannot be carried 
as far analytically. 

In general the pairing interaction V can have off-diagonal as well as diagonal terms in 
the L basis, and in general all channels will be coupled. If the interaction has the symmetry 
of the Fermi surface then the integral equations only couple different channels L and V 
which transform according to the same irreducible representation. However, the subsequent 
matrix can be diagonalized with respect to the indices of the same representation resulting 
in a new set of basis functions which are linear combinations of the old basis functions of 



the same representation in different L channels Thus we can then write the interaction 
as a diagonal matrix in the new basis functions which still has a general structure for each 
representation within each new channel. Thus we can write V^'y = Vi^^L,v^m,m' ■ This 
allows us to reduce the infinite series of coupled integral equations to a limited subset that 
can be handled analytically. The elements of the expansion will be dominated by a single V™ 
component corresponding to the L pairing channel symmetry m, and the other components 
represent admixtures of channels (the smaller eigenvalues in the gap equation) with different 
pairing symmetry. If the gap representation is one dimensional (all representations of the 
jjAh g roU p except Eg, which is two dimensional), then all the other V^'s are set to zero 
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and there will only be collective modes connected with the broken U(l) gauge invariance. 
Otherwise other collective modes can be present as well. 

We now define new vertices T^(iu) = 7™ + V™8™(iu) and expand 6™(iu) in spin quater- 
nions, 

3 



[IUJ] 



(C8) 



i=0 



^=2 



(k) 



The index m stands for the representations of the D 4h point group, and are defined through 
the basis functions as follows 

B lg (T+) (m=l), 
A ls (r+) (m = 2), 

B 2g {Tt) (m = 3), (C9) 
E g (Tt) (m = 4), 
E g (Ft) (m = 5), 

where the Tf corresponds to the notation of Sigrist and Rice [p5f] . It can be shown that 
the coefficients m of Pauli matrices tq^ satisfy homogeneous equations and thus vanish 
while the remaining coefficients satisfy the following coupled integral equations, 



( 75 (^s k y)> 
h h 

k x k Z) 

h h 
rvyrv z , 



5i 2)m M = £ {V^) m (iu)C^ m '(iu) + i(7^' + V^5f} m \iu))A^\iu)}, (CIO) 



L':m' 



4 3)m (^) = E {iV$ 5$™' [iu)A™$ {iu) - ( T £ + V^) m (iu))C^r'(iu)}, 

L>;m> 



where the functions A, C are given by 



A L,L> (M 

C L . L > («*>) 



^(k)^(k)A(k,ia;)>, 



The spectral functions are defined as 

A(k,iu) = -iA(k)uN F \ d£ 

C + (k,iu) = —Np f d£ 



1 l-2f(E) 
AE 2 iu-2E 
A(k) 2 l-2/(£) 



[IUJ 



-IUJ) 



2E 2 iu - 2E 



[iuj 



-iu), 



-IUJ). 



(Oil) 



(C12) 



(C13) 



Here / is a Fermi function and (iu — > — iu) denotes additional terms which differ only in the 
sign of iu. Analytically continuing to the real axis by letting iu — > uj + iO, the £ integration 
can be performed analytically and we obtain 
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A(k,u) = A(k) F(k,u), 
C +(k,uj) = ^^-F(k,uj) 
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(C14) 



C-(k,u) 



-^ + -F(k,u) 
V(k) 2 y 1 ' 



with 



N F 



F(k,u) 



VA(k) 2 -(w/2)2 
Np 

2 v /(^/2) 2 -A(k)2 



arctan 



2v/A(k)2-(o;/2)2 
Z7T + log 



for A(k) 2 > 



^/2-y/(a;/2)2-A(k)2 



and 



.^/2+v/(a;/2) 2 -A(k) 2 



:t) 2 , 

'aA2 



for A(k) 2 < (|) , 



(C15) 



(C16) 



V(k) " * Jo 2E(k) 

given by the BCS gap equation. The function F is closely related to the Tsuneto function, 
Eq. (11). 

The integral equations are still general for a charge density-like vertex and the symmetry 
of the interaction and gap remain undefined. We now restrict our attention to the case of 
o?-wave interactions such that only V^ =2 7^ an d other terms corresponding to interactions 
in higher angular momentum channels are discarded. Dropping the L = 2 subscript by 
denoting V£L 2 by V m and 5^=2™ by £m' 3 \ the integral equations simplify to 



<J< 2) < 



ILO) 



(C17) 



^ 3) ( 



lOJ) 



hm,m 
'L=2,L=2 



(to)+id^(to)A^ 



=2,L=2' 



(C18) 



These equations are still general to any <i-wave pair state. 

We now specifically work with a B lg gap, A(k) = A (/c 2 — k^) (the Tg representation 
|25[), noting that similar conclusion can be drawn for other choices of energy gaps within 
the L = 2 subgroup of the -D 4h point group. The coupled integral equations represent 
10 equations for the 10 unknowns S^,i G {2,3}, and the solution can be obtained by 
diagonalizing a 10 x 10 matrix. In aiding to solve the coupled integral equations, it is useful 
to examine the selection rules of the spectral functions A, C ± for various subgroup indices 
within the L = 0, 2 channels. In particular we note that 
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s-i±m,m' r ^±4,4 _ x-y±5,5 syiO,m r 

°L,L ~ O m>m /, <~'L=2,L=2 — °L=2,L=2' <^L=0,L=2 ~ °m,2, 



A lS,L=2 ~ ^m,m', for 771, 77l' ^ {2, 3}, -4l=2,L=2 ~ <W, ^L=2,L=2 ~ 5 



m.l ; 



^L=0,L=0 — 0' ^L=2,L=0 ~ ^m,l) ^L=2,L=2 — ^L=2,L=2- (C19) 

These theorems hold for the case of a B lg gap only. 

Handling the full gauge invariant response amounts to solving a matrix equation to 
identify the collective modes in each channel. Solving these equations allows one to identify 
the position of the collective modes by locating the zeroes of the real part of the denominator. 
We first simplify our notation by defining 

L,m' 

/£M = - E il cIXl^)- (C20) 

L,m' 

We find that the m = l(B lg ) and m = 2(Ai g ) channels are coupled due to the fact that a 
d x 2_ y 2 energy gap squared has a component which has a finite overlap with the A lg channel 
which is isotropic within the x — y plane. Solving the integral equations we obtain for the 
m — 1, {Big) and m — 2, (A lg ) channels, 

X (2 )r , _ J 1 A (iuj) + V 2 Al' 2 2 (iuj)5 ( i\tuj) 
°l { luj )—l - n , 

1 - 14C 2)2 

r(3),. v fM 

o{ yioj) = - - - - 



1 i vr+^dtA i ^i^jMi 2 ' 

l + l/iG 2)2 + ^c-vi^ 



(2) Jj{tuj) + V l A l ^{tu J )5f\iu J ) 

1 - 1/262,2 (ioj) 

di' «w = , 10, — ■ 621 

1 + 1/262,2 + IZ^c^M 

For the m = 3(S 2g ) and m = 4, 5 (two channels, we find no coupling and obtain 



fCl 



4 3) M = fjL Mr, r , (C22) 



jf(2)/. v J 4 A 5 (^) + ^4,5^2;2 5 ' 5 (^)^4 3 5(^) 

5 4M^) = « " T/ ^-4,4;5,5/- \ 

1 - V 4 , 5 6\2 (ttj) 
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where V^, 5^, C ±4 ' 4;5,5 , and A 4,4 ' 5,5 stand for the interactions, renormalized vertices, and 
spectral functions in the m = 4, 5 channels, respectively. These equations represent the full 
channel dependent renormalization of the Raman vertex due to d— wave pairing interactions. 

We now are in a position to identify the collective modes in the 10 channels (5 real and 5 
imaginary) by locating the zeroes of the denominator in each channel. While the functions 
f A ' c remain general and are determined by the magnitude of the Raman vertex they do 
not enter into the denominators of the expressions for the renormalized vertices and thus 
are irrelevant for the position and broadening of the collective modes. The modes located 
by the zeroes of the denominator are of two types, namely, massive and massless. The 
massless (Goldstone) modes are a consequence of the spontaneously broken U(l) continuous 
symmetry by the interactions, while the massive (optical) modes are generated as well. Using 
Eqs. (C21-C23), in each channel we locate the zeroes of the real part of the denominator in 
each channel to find the position of the collective mode and evaluate the imaginary part of 
the denominator at the position of the collective mode to determine its broadening. 

Our results are summarized in Table 1. We see that the Goldstone mode (u c = 0) 
appears in the A\ g channel and the other modes are all massive. The B\ g modes lie beneath 
the maximum energy gap while the others appear above 2A . These excitonic modes are 
damped considerably due to the existence of quasiparticles from the presence of gap nodes 
which provide decay channels to damp the massive modes. This is in contrast to excitonic 
modes in s-wave superconductors pl| , p3| . However, we do observe in Table 1 that the mode 
position rapidly decreases to lower frequencies with depreciating residue as the coupling 
stregnth V m /Vg lg is reduced, which is similar to the BCS case. Below a critical coupling 
strength V m /VB lg ~ 0.8 the collective modes disappear altogether and thus have little impact 
on the low frequency behavior of the Raman correlation function for small couplings. 

We now reconstruct the Raman response and determine which of the collective modes 
couple to the Raman vertex. Putting in Eqs. (C21)-(C23) into (CI), we can express the full 
gauge invariant Raman response function as 

X£(q=0,ia;)= (C24) 
E(E7^?'' m 'M + Vm' 6$ (w) Ct™' m ' (iu) - iV m ,8%{iu)A™f{iu>)). 

m> L> 

Taking only the L = 0, 2 terms of the Raman vertex into account (see Section III), we carry 
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out the summation over V and ml to obtain the Raman spectrum in each channel, and then 
take into account the long-range Coulomb screening via [^l[ 



XsM - o, iu>) = x 7 , 7 M - X7 ' 7 ^°7\ > (C25) 

where Xa,b denotes Eq. (CI) with the vertices A(k) and B(k), respectively. 

We find that for the density channel (L = 0), x(q = 0, u) — 0, which is a restatement of 
particle number conservation and bears out the gauge invariant nature of the theory. This 
simply restates that intercell charge fluctuations couple to the long-range Coulomb forces to 
be completely screened for q — > 0. After lengthy but trivial algebra we obtain the compact 
results for the intracell fluctuation contributions. For the B\ g channel we obtain, 

l + ViC Blg {iuJ) 



and for the B 2g channel, 



XB 2fl M = 2( 7 f) 2 — , (C27) 



and for the two E g channels, 



^-*«^+ifisp (C28) 



with the functions 

C Blg (iu) = C$\iu) + V 2 '^^j 2 , (C29) 

1 - V 2 C 22 ' (zu) 

C Bag (iw) = C$> 3 (iu), (C30) 

I ^4(5),4(5) , • n ,2 

c,„M = cjf >«M + ^-^-S-. (C3i) 

1 - ^4(5)02,2 

These channels are unaffected by screening since the intracell fluctuations lead to no net 
charge transfer for these symmetries. The expression for the fully symmetric A\ g channel 
which contains both inter-and intracell fluctuations expression is much more complicated 
due to Coulomb screening, 



I — 32. y (J® (iij) 
1 + V 2 C Alg (iuj) 
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Here the functions CU lg , C° Mn are defined as 



Ao'o(iuj) 



C Alg (iu) = Ctfiiuj) + V 1 1 "%Zr r , (C33) 

1 - ViC 2 2 ' [iu) 



and 



CU lfl M = C 2 ,o M + V. n • (C34) 

Examing the structure of the denominator of Eqs. (C26)-(C31), we see that the massless 
gauge mode (Goldstone mode) in the A\ g channel drops out of the Raman response, and 
only provides screening via Eq. (C25). The massive modes do appear in the B\ g and E g 
channels while the massive mode found in the B 2g channel does not couple to a Raman 
probe. Linearizing the denominator around the position of the collective mode we calculate 
the residue Z of the mode to be Z Blg = uj log(A /o; c ) and Z Eg = u c , where u c is the 
position of the collective mode in each channel. Therefore, putting all our results together 
for the collective modes, we can argue that the collective modes are of little relevance to 
electronic Raman scattering due to the fact that l)the modes only exist at beyond a large 
strong couplings threshold, 2) the residue of the modes decreases the lower the position of 
the collective mode, and 3) that the broadening of the collective modes are substantial. 

However, of course the final-state interactions themselves can affect the spectrum |2T,23|. 
Therefore we display our results in Figs. (8a-c) for the entire spectrum in the B lg , B 2g , 
and A lg channels for different ratios of the parameter V m /VB 1 , where Ve l9 is the pairing 
interaction of d x 2_ y 2 symmetry. The response for the E g channels are similar to the B 2g case. 
We note that the interactions have only a minor affect on the spectra in the B\ g channel, 
only changing the cusp behavior near 2Ao, (this demonstrates how the finite ^-dispersion 
cuts off the logarithmic singularity in Fig. 1 for the 2-D Fermi surface), while leaving the 
peak position and low frequency behavior unchanged. The interactions have more an effect 
in the B 2g and A\ g channels due to the fact that the peak of the spectra are very broad. 
Therefore the interactions shift the peak position upwards in frequency along the top of the 
broad hump of the spectrum, from 1.3 to 1.5Ao for the B 2g channel, and from 0.6 to O.8A0 
in the A\ g channel. Again we note that the low frequency behavior remains unchanged. A 
similar effect is seen in the E g channels. In particular, the massive mode in this channel is 
entirely damped leading to no drastic changes in the spectrum. 

We close this Appendix by addressing the collective modes for a superconductor with 
a cylindrical (2-D) Fermi surface. In a 2-D system the only interactions that appear at 
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the L = 2 level are the B\ g and B 2g channels, while the A\ g and A 2g channels appear at 
L — 4. Similarly, there are no E g channels for a system without dispersion in the z-direction. 
Therefore, if we are only concerned with o?~wave interactions and not those of higher angular 
momentum channels, then we can set the matrix elements Va 1b = Ve (1,2) = in Eqs. (C21- 
C23) and the equations simplify greatly. Since then there will be no channel mixing of the 
A\ g channels into the response, then therefore there can be no collective mode since no 
term appears in the denominators which has the form 1 — V m /VB 1 which arises through the 
function C~ (see Eqs. (C21-C23)). Therefore we can conclude that collective modes can 
only appear when the interactions are included in higher order L channels for a 2-D system. 
Furthermore, the vertex corrections produce minor changes in the spectra only for the B\ g 
and B 2g channels. The changes are similar to the changes shown for the response functions 
evaluated on a spherical Fermi surface (apart from cutting off the logarithmic divergence of 
the Bi g response at the gap edge), and thus are not of major importance. We can therefore 
neglect the collective modes entirely and the effect of vertex corrections and simply use the 
bare bubble for the Raman response. This completes the purpose of the Appendix. 
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FIGURES 

FIG. 1. Electronic Raman response functions evaluated for d x 2_ y 2 pairing on a cylindrical Fermi 
surface for various symmetries as indicated. All vertices have been set equal to 1. 

FIG. 2. Ratio of the low frequency Raman response in the superconductor with d x 2_ y 2 pairing 
to the normal metal for the temperatures indicated. Note the slow decrease of the A\ g and B<i g 
channels with temperature. 

FIG. 3. Response functions for various symmetries for s + id pairing with A s /A^ = 0.25. All 
vertices are set to 1. 

FIG. 4. Response functions for various symmetries for anisotropic s-wave pairing with 
Aq/Ai = 0.25. All vertices are set to 1. 



FIG. 5. Comparison of the theory to the experimental data taken on BSCCO from Ref. [12] 
using d x 2_ y 2 pairing. The parameters used are defined in the text. 

FIG. 6. Comparison of the theory to the experimental data taken on YBCO using d x 2_ y 2 
pairing from Ref. [p!H| . The parameters used are defined in the text. 

FIG. 7. Comparison of the theory to the experimental data taken on single layer Thallium 
cuprate using d x 2_ y 2 pairing from Ref. The parameters used are defined in the text. 

FIG. 8. Effect of vertex corrections on the Raman response evaluated for d x 2_ y 2 pairing on a 
spherical Fermi surface for (a) B± g , (b) E>2 g , and (c) Ai g channels (using 70/72 = 2, other vertices 
set equal to 1). The E g results look identical to the B2 g spectra. Here we have used Vb 13 = 0.2, 
and the values of V m /VB lg are indicated in the upper right hand corner of the Figure. 



42 




Devereaux and Einzel, Fig. 1. 
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TABLE 1 

Position cj c /2A and broadening r c /2A of the pole in for each channel. 
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Goldstone 


Goldstone 
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1.09 


r c 
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0.20 
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